Cluster persistence in one-dimensional diffusion— limited cluster-cluster aggregation 
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The persistence probability, Pc(t), of a cluster to remain unaggregated is studied in cluster- 
cluster aggregation, when the diffusion coefficient of a cluster depends on its size s as D(s) ~ s 1 . 
In the mean-field the problem maps to the survival of three annihilating random walkers with time- 
dependent noise correlations. For 7 > the motion of persistent clusters becomes asymptotically 
irrelevant and the mean-field theory provides a correct description. For 7 < the spatial fluctuations 
remain relevant and the persistence probability is overestimated by the random walk theory. The 
decay of persistence determines the small size tail of the cluster size distribution. For < 7 < 2 the 
distribution is flat and, surprisingly, independent of 7. 
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I. INTRODUCTION 

Aggregation models are useful in describing various 
phenomena from chemical engineering, material sciences, 
atmosphere research to even astrophysics |Q, ^ ||. One 
general property of these models is that they lead to dy- 
namic scale-invariance: when all the lengths are scaled 
by the characteristic length, the system looks the same 
at different times. Lately, studies of first passage prob- 
lems (J] under the name persistence § |, § have shown 
that not necessarily all the properties of a dynamical! 
scaling system are characterized by a single scale 
Here we address the probability of a cluster to remain 
intact in an aggregation system and show how this quan- 
tity and the associated length scale relate to the physi- 
cally relevant issue of the shape of the cluster size distri- 
bution. 

In an aggregation system one can define many first- 
passage problems and related quantities ||. We study 
the probability that a cluster has not aggregated with 
any other one before time t jl0| . This probability is called 
cluster persistence and denoted by Pc(t). Similar prob- 
lems considering uninfected walkers in one-dimensional 
reaction-diffusion systems |ll| and Potts model jl2| have 
recently been shown to display interesting behavior. We 
concentrate on diffusion-limited cluster-cluster aggre- 
gation (DLCA) in one dimension, where the dynam- 
ics is dominated by spatial fluctuations For high 
dimensional systems these may be neglected, and on 
the mean- field level, valid for dimensions higher than 
the upper critical dimension, aggregation is well under- 
stood 

The DLCA model is defined so that the nearest neigh- 
bor occupied sites in a lattice are identified as a cluster. 
Each cluster diffuses with a size dependent diffusion con- 
stant, D(s) ~ s 7 , where 7 is the diffusion exponent. If 
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a cluster collides with another one, the two clusters are 
irreversibly merged together and the aggregate diffuses 
either faster (7 > 0) or slower (7 < 0) than either of the 
colliding clusters. While the size independent diffusion 
(7 = 0) is exactly solvable in one dimension, it forms a 
marginal case between two completely different aggrega- 
tion mechanisms |l7j . We study here the more physically 
interesting problem with 7 7^ 0. 

The aim is to study the dependence of cluster persis- 
tence on the diffusion exponent 7 and extend the study 
presented in Ref. |Io| . We also pay attention to the 
random-walk (RW) problems that ensue as on a mean- 
field level the problem is reduced to the survival of three 
annihilating random walkers. While the 7 = case is 
readily solvable by various methods [ll| , already the 
case of three annihilating particles with unequal diffusion 
constants is rather involved p(| . Here 7^0 leads to time 
dependent diffusion coefficients, and we derive a Fokker- 
Planck (FP) equation for the survival of these particles. 
For 7 > its analysis yields an algebraically decaying 
survival probability, P sn i V (t) ~ t^ 9 ™^. The survival 
exponent #kw is discontinuous and non-monotonic as it 
is given by 9rw(i) = 2/(2 — 7) for < 7 < 2 and 
$rw(0) = 3/2. The numerical comparison of the sur- 
vival and persistence probabilities validates the theory 
and hence Pc(t) ~ t~ e ° with 9c — 9rw- 

For 7 < simulations show that both the survival 
and persistence probabilities decay stretched exponen- 
tially as exp(— Ct"). The Fokker-Planck equation is not 
amenable to analytic analysis, so we use a Lifshitz tail 
argument to understand the survival. Such heuristic 
arguments and numerics suggest a stretching exponent 
Pyw(i) — ^7/(4— 27). The Lifshitz tail argument indi- 
cates that the exponent is affected by the fluctuations in 
the motion of the particles that neighbor the surviving 
one. These are taken into account only approximately in 
the mean-field theory and for the DLCA numerics gives 
Pc = —27/(6 — 37). A closer examination reveals that 
also the distance distribution between the particles sur- 
rounding a surviving one in the mean-field model scales 
in a different way than the corresponding distribution of 
the DLCA. 
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In addition, we show how the cluster persistence is 
related to the cluster size distribution. To clarify the 
connection, consider the dynamic scaling in DLCA. Both 
simulations and experiments show that the cluster size 
distribution n s (t) (the number of cluster of size s per 
lattice site at time t) scales as |lj 



,(t) - S(ty 2 .f 



S(t) 



(1) 



where S(t) ~ t z is the average cluster size and the scal- 
ing limit, s — > oo and S(t) — > oo with s/S(t) fixed, 
is taken. In one dimension the dynamic exponent z = 
1/(2 — 7) [Hi HI' ^ or 7 — the cluster size distribution 
is broad in the sense that the scaling function behaves 
as f(x) ~ x~ T as x = s/S(t) — ► 0. For 7 < the 
scaling function is bell-shaped and f(x) ~ exp(— AaH") 
for x — > 0, where ^4 is a constant. To determine the 
polydispersity exponent, r, which characterizes the num- 
ber of small clusters, is non-trivial even on a mean-field 
level [[l4| ^3) whereas the similar exponent n readily fol- 
lows from scaling analysis p5| . All the exponents z, t, 
and fj, are universal, i.e., they do not depend on the fine 
details of the model. They can, and it is natural to expect 
that they do, depend on the diffusion exponent 7. 

One of the main results of this paper is that the expo- 
nents describing the decay of the cluster persistence are 
related to these universal exponents as 



"c 

Pc 



(2-t)z 
\li\z. 



(2a) 
(2b) 



Quite unexpectedly, the polydispersity exponent is a con- 
stant, t = 0, for < 7 < 2, but discontinuous since 
t(7 = 0) = —1. The reasoning leading to the rela- 
tions (|2a|) and (2b) is universally applicable, so that the 
behavior of the tail of cluster size distribution might be 
tackled through cluster persistence in other models, too. 

The outline of the paper is as follows. In section || the 
mean-field random walk theory is formulated and the as- 
sociated Fokker-Planck equation is derived. Section |l| 
starts by describing the simulation methods. Thereafter 
the mean-field theory is validated for 7 > by compar- 
ing the survival probability obtained from the analysis of 
the Fokker-Planck equation to the simulation results of 
both the random walk system and the DLCA one. For 
7 < a similar comparison shows the effect of spatial 
fluctuations, and the stretched exponential decay of the 
survival probability is explained using a Lifshitz tail ar- 
gument. Section IV concentrates on the relation between 
the persistence and the small size tail of the cluster size 
distribution. The paper ends with conclusions in sec- 
tion |v|. 



persistent one). The cluster in the middle will be per- 
sistent until it collides with one of the neighbors. Af- 
ter this the two remaining clusters would contribute to 
persistence only by increasing the mass of the clusters 
surrounding another persistent cluster. This is negligible 
at late times, since the persistent clusters will be sep- 
arated by many non-persistent ones, i.e. t 9c > t z . In 
other words, the correlations in the system grow only as 
t z and each persistent cluster is asymptotically indepen- 
dent. Thus it is sufficient to consider only one persistent 
cluster and its two neighbors. 

The collisions of the surrounding clusters make them 
bigger and increase or decrease the diffusivity. We make 
the mean-field approximation that each cluster neigh- 
boring a persistent one will grow as an average cluster 
does. Hence, we replace the true process, where the sur- 
rounding clusters collide at some discrete times U, by a 
continuous one, where the surrounding clusters grow as 
S(t). As -D(s) ~ s 7 these clusters will diffuse with time- 
dependent diffusion coefficients. In the following analysis 
we will ignore the possible early time crossover effects in 
the growth of the average cluster size and the diffusion 
coefficients of the clusters surrounding a persistent one 
are taken to follow a true power-law at all times. This 
will only affect the early time behavior. 

The finite extent of clusters is irrelevant for cluster per- 
sistence and we will consider the three clusters as point- 
like particles from now on. Let Xi(t) (i = 1,2,3) de- 
note the positions of the particles at time t such that 
x\(0) < x 2 (0) < 2:3(0). The motion of these is described 
by the Langevin equations 



(3) 



with Gaussian white noises (£i(t)) = and (£,i(t)^j(t')) = 
2T>i(t)Sij5(t — if). The overdot denotes derivative with 
respect to time and the brackets an ensemble average over 
different realizations. The diffusion coefficients of the 
particles read as T>i(t) = T> 3 (t) = D\V Z and T> 2 (t) = D 2 . 
The meaning of a time-dependent diffusion coefficient, 
say T>i (t), is simply that the particle 1 will follow a simple 
diffusive motion with a diffusion constant D\ in the time 
scale 



Ti(t) 



&'V 1 {t')/D 1 =f< z + 1 /{ 1 z + l). (4) 



As we are interested in the survival of the middle par- 
ticle (x 2 ), the process terminates when cither x\ = x 2 
or x 2 = X3. It is convenient to consider the distances 
between the particles: xi 2 (t) = x 2 (t) — Xi(t) > and 
^23 (i) = %3{t) — x 2 (t) > 0. These obey similar Langevin 
equations 



II. MEAN-FIELD: REDUCTION TO A THREE 
PARTICLE PROBLEM 

The two clusters surrounding a persistent one will grow 
when they collide with other clusters (but not with the 



iia(*) = r u (t) 

±23 (t) = r 23 (t), 



(5) 



where (T 12 (t)) = {T 23 (t)) = and (r 12 (t)r 12 (i')) = 
(r 2 3(«)r 2 3(£')) = 2(L>2 + D^^Sit - t'). The two noises 
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are correlated as the motion of the middle particle af- 
fects both distances: (T 12 (t)r 23 (t')) = -(6(*)6(*0) = 
— 2D 2 8(t — t') 7^ 0. For 7 > the noise correlations be- 
come asymptotically irrelevant, which is not the case for 
7 < 0. 

To proceed, we transform Eqs. ((H) to a Fokker-Planck 
equation for the probability density p(cci2, X23] t) of the 
two distances at time t. Due to the mutual correlations 
this is easiest to do by computing the drift and diffusion 
coefficients from their definitions 

A = A lim o — (xi (t + At) - Xl 0)) 

and insert these to the general Fokker-Planck equa- 
tion M 



dp 
dt 
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d 2 



dxidx-j 



(6) 



A straightforward calculation gives 

m- {D2 + D ^ ) {dxT 2 + dxf 3 



ID, 



d 2 P 



'dxi 2 dx 2 3 



(7) 



The initial condition is now p{x\ 2 , X23] 0) = S(x\ 2 — 
x i 2 )fi( x 23— x 23 ), where x\ 2 = £12(0) and = £23(0) are 
the initial distances between particles. The termination 
of the process when two particles collide gives absorbing 
boundary conditions along the axis, i.e., p(x\ 2 ,Q;t) = 
and p(0, x 2 3\ t) = for all times t. 

Thus the original many body problem has been re- 
duced to the survival of three annihilating random walk- 
ers. Given that one can solve Eq. (Q) with the appropriate 
boundary conditions, the survival probability of the mid- 
dle particle (which corresponds to the persistent cluster) 
can be obtained as 



between neighboring clusters or random, in which case 
each site is independently filled with probability <j). The 
persistence exponent is independent of the initial dis- 
tribution but the early time behavior of the persistence 
probability depends on it ||. 

In the dynamical evolution a cluster is selected ran- 
domly and time is increased by l/[N(t)D max (t)]. Here 
N(t) denotes the number of clusters and D max (t) is the 
maximum of the diffusion coefficients of all the clusters 
at time t. The cluster is moved one lattice spacing with 
cluster size dependent probability D(s)/D max (i). If the 
cluster collides with another one, the clusters are irre- 
versibly aggregated together and the values of N and 
-Dmax are updated. Then a new cluster is selected and 
the above procedure is repeated. 

The three particle simulation is similar to that of the 
DLCA. Initially the distance between particles is Zo- A 
particle i £ {1, 2, 3} is selected randomly and it is moved 
a distance a either to the left or to the right with proba- 
bility Vi{t)/V>{t). Here £>>(£) = max{I?i(t), D 2 ,V 3 (t)} 
is the maximum of the diffusion coefficients of the three 
particles at that time. The distance a is set to corre- 
spond the lattice constant of the DLCA simulations, i.e., 
a = 1 . Irrespective of the movement of the particle time 
is increased by l/[3X>>(t)] and the time-dependent dif- 
fusion coefficients 2?i(i) and T> 2 (t) are updated to new 
values. This procedure continues until a collision occurs. 
Figure [l] shows examples of configurations that survive 
for a long while for negative and positive values of the 
diffusion exponent. 

The faster the survival probability decays the more 
computation time is used in simulating systems, which 
terminate at early times. In order to sample efficiently 
the long living, interesting configurations we use a cloning 



method [|25|, |26| for the three particle simulations when 
7 < 0: At times tj we make rij copies of all the systems, 
which have survived upto this time. Typically simula- 
tions are averaged over 5 x 10 6 initializations and the 
system is copied at times 10, 10 2 , and 10 3 with 2 x 10 3 , 
5 x 10 4 , and 10 4 copies, respectively. This enables us to 



-Psurv(i) = / dxi2 / dx 2 3 p(xi 2 ,x 2 3;t). (8) reach probabilities less than 10 
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When the survival probability decays algebraically, 
isurv(i) ~ t~ e ™ , the associated exponent 6>rw is called 
the survival exponent. 



B. Size-Independent Diffusion (7 = 0) and 
Crossover Behavior 



III. COMPARISON OF THE SIMULATIONS 
AND THEORY 

A. Details of Simulations 

The DLCA simulations are done on a lattice of size L 
with periodic boundary conditions. Concentration cf> of 
sites is filled with particles and nearest neighbor particles 
belong to the same cluster. The initial distribution is ei- 
ther monodisperse, n s (0) = <5i. s , with equal distances, Iq, 



When the diffusion constant of a cluster does not de- 
pend on its size, i.e. 7 = 0, an exact solution is possible 
as the collisions of the clusters surrounding a persistent 
one with other clusters do not matter ^7j. For the same 
reason the mean-field approximation becomes exact and 
reduces to an old problem of the survival probability of 
three similar annihilating random walkers 1 0| . The per- 
sistence and survival exponents attain the value 3/2. 

This result can also be obtained from the equation d?]) 
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FIG. 1: Visualization of the three particle system when the initial distance lo = 10 for a) 7 = —2 [z = 1/4) and b) 7 = 1 (z = 1). 
The probabilities of these configurations are of order 10 -8 and 10 -4 , respectively. At the final time t he rat io D2/(Dit~ lz ) is 
about 10 2 in (a) and 10 -3 in (b). The dashed lines in (a) show ^"-behavior with a = 0.375 (see Section [II D for details). 



which for this particular case simplifies to 



C. Validation of the Mean-Field Theory (7 > 0) 



dp 
dt 



d 2 p d 2 p 



y dx\ 2 



dx\ z 



d 2 P 



8x128x23 



(9) 



where we have taken D\ = D2 = 1/2. A coordinate 
transformation x = [xn + £23)) 2/ = { x \2 — a^23)/v / 3 re- 
duces this to a diffusion equation 



dp 
dt 



dx 2 



d'p 
dy 2 



(10) 



with the boundary condition p = along lines y = 
±a;/v3. This corresponds to a two-dimensional wedge 
of angle = 7r/3, in which the survival probability de- 
cays as t~^ 2e ~ t~ 3 / 2 g. 

It is also interesting to know how the asymptotic 
regime, where P surv (t) ~ i~ eRW , is reached. In the case 
D-y = D% = D3 = D (7 = 0) with the initial distances 



between particles being x\ 2 



1 23 



Iq the solution in- 



cluding the first correction to scaling is given by |19L [27 



r(t) 



4V2T \Dt) 



3/2 



72 
'0 



3 

16 Dt 



(11) 



The correction becomes negligible for times much larger 
than the crossover time t CT = 3Zg/(16D). For 7 ^ the 
corrections go in powers of the ratio of the diffusion co- 
efficients, D2I D\V Z . For 7 > this is demonstrated in 
Appendix ^ and for the corresponding two particle prob- 
lem it may be shown exactly (see Appendix |b|). There- 
fore the crossover time depends on 7 as t cr ~ 7~' 2 tJ/ItI a 
where the constant r ~ 30 according to simulations. As 
t CI diverges for I7I — * 0, we can expect that the asymp- 
totic scaling regime can be reached in simulations only 
for relatively large values of |7|. 



We have not been able to solve equation (0) exactly. 
The reason is that the absorbing boundary conditions to- 
gether with the two time scales appearing in the problem 
make the standard methods (Laplace or Fourier trans- 
forms; polar coordinates) unapplicable. Nor is it pos- 
sible to transform the equation to a diffusion equation 
with simple enough boundary conditions. However, the 
full solution is not needed for the determination of the 
survival exponent since this is given by the leading large 
time behavior when t — > 00. It would only provide us in- 
formation about the crossover effects, which according to 
our analysis (see Appendix ^) and the numerical simula- 
tions (see below) are rather pronounced when 7 is close 
to zero. 

A change of variables x = (xi% + X23)/V%, y = (x\2 — 
a ; 23)/v / 2 transforms Eq. (0) to 



dp 

dt 



cfp_ 

dx 2 



[D X V Z 



(12) 



with the boundary condition p = along y = ±x, i.e., 
a wedge of angle 6 = ir/2. When 7 > the con- 
stant term is negligible at long times {D\t lz ^S> D 2 ) 
and the diffusion becomes isotropic. This can be shown 
by directly solving equation (^) and analyzing the large 
time behavior of the solution (Appendix [X]) . A change 
to the time scale T\ [see Eq. (0)] transforms Eq. ( |l2] ) 
to the form of Eq. (10) and the survival probability 
P SU rv(<) ~ T-" /2& ~ Tf 1 ~ t-^ z \ As z = 1/(2 - 7 ) 
the survival exponent = 2/(2 — 7) = 2z. 

The approximation of neglecting the constant term in 
Eq. (^) corresponds to a complete separation of the time 
scales, i.e., to a situation, where the middle particle is 
at rest {D2 = 0). Thus for 7 > one could simply 
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determine the survival exponent by considering two in- 
dependent random walkers with a fixed absorbing bound- 
ary in between [compare to Fig. |l| (b)]. In other words, 
the motion of the "slow" particle becomes asymptotically 
irrelevant. This can be exactly shown for the the corre- 
sponding two particle problem (Appendix |b|). 

Figure ^ compares the survival and persistence prob- 
abilities. The initial distances between particles in the 
random walk simulations are set to be the same as in 
the DLCA. The probabilities decay algebraically at large 
times and the only difference in the decay is between the 
amplitudes. This is to be expected as the transient ef- 
fects of the growth of the average cluster size are not 
taken into account in the random walk picture. 

The inset shows local exponents, i.e. logarithmic 
derivatives of the probabilities, which converge to the 
value obtained from the Fokker-Planck equation, (9rw = 
2/(2-7) for 7 > and 9 RW = 3/2 for 7 = 0. The asymp- 
totic regime is reached only for 7 = and 7 > 0.5. In the 
latter region the local exponents saturate, when the ratio 
of the diffusion coefficients is of about 30. For example, 
for 7 = 0.25 this would corresponds to t cr 2 x 10 10 
which is beyond the time reached in simulations. 

Note that the persistence exponent is discontinuous 
and nonmonotonic at 7 = 0, i.e., 3/2 = 9c(0) > 
$c(0 + ) = 1. This seems first counterintuitive since mak- 
ing some of the clusters to diffuse faster helps others to 
survive longer! On the other hand, as time elapses a 
persistent cluster becomes slower as compared to an av- 
erage one. In this way it eventually adopts the optimal 
strategy S| by becoming stationary. 



D. Fluctuation dominated persistence (7 < 0) 

For 7 < the diffusion of the clusters surrounding a 
persistent one slows down. Consider the random walk 
picture and proceeding similarly as for 7 > above. 
Fixing now particles 1 and 3 would lead to an inter- 
val of fixed length and hence to an exponentially de- 
caying survival probability. However, simulations show 
that the survival decays stretched exponentially in time, 
P S uTv(t) ~ exp(— Crw^ ,3rw )- Furthermore, as will be 
shown below, although the surrounding particles become 
slower, their motion can not be neglected even at the long 
time limit. This is a collective effect and in clear contrast 
to the exactly solvable two particle case, where the fast 
particle eventually dominates the survival (Appendix^). 

In figure || we plot — ln(P surv (t)) vs. t on a log- log 
scale so that a stretched exponential decay corresponds 
to a straight line with a slope /3rw The final slope is 
independent of the initial distance between particles, and 
thus the stretching exponent is universal. 

Figure || shows the location distribution p(x 3 ; t) of the 
particle 3 (the one for the particle 1 would be the same) . 



It scales as 



p(x 3 ;t) = t~ 



bt a 



(13) 



implying that although the distribution widens as t z , the 
expectation value of the distance from the origin grows 
as bt a with a non-trivial exponent z < a < 1/2 [see 
Fig. (a)]. The scaling is similar to the reaction front 
in the originally separated reaction-diffusion system A + 
B — > C, where the reaction zone becomes sharp at late 
times, i.e. z < a 29, |3(|. It is striking that the scaling 



function g(y) is within the numerical accuracy a simple 
Gaussian. 

The consequence of Eq. ([[3]) is that the average dis- 
tance between the particles 1 and 3 grows [see Fig. |l| (a)]. 
If it would grow deterministically as t a , with a < 
1/2, the survival probability would decay asymptoti- 
cally stretched exponentially with the exponent /3dct = 
1 — 2a |3^[. For example, for 7 = —2.0 the numerics gives 
a rough estimate a 0.36 and 1 — 2a ~ 0.28, which is 
in reasonable agreement with the numerically obtained 
stretching esponent /3rw ~ 0.25 (see inset of Fig. ||). 

To understand the origin of the new length scale t a 
the next logical step is to try to take the length fluctua- 
tions of the interval into account. We make this using a 
Lifshitz tail approach W. It is based on the assumption 
that the main contribution to the survival is provided by 
extreme configurations, where the particles surrounding 
the surviving one have diffused far apart from each other. 
We write the survival probability as 



-^surv (^) 



dl P(l;t)Q(t\l), 



(14) 



where P(l;t) is the probability distribution of the inter- 
val lengths I = X3 — x\ around a surviving particle at 
time t and Q(t\l) ~ P 1 exp(— n 2 Dt/l 2 ) is the survival 
probability of a particle in an interval of length I ||, [52| . 
In order to make progress, we need to know the large I 
behavior of P(l; t). It scales similarly as p(xz; t) 



P(l;t)=t~ z G 



I - 2bt a 



(15) 



where the large y tail of G(y) is Gaussian as the position 
distributions of particles 1 and 3 are Gaussian. Although 
it is irrelevant in what follows, the small y part of G(y) 
decays faster than the large y tail due to the restriction 
x 3 > Xx. 

Denote the variance of the Gaussian tail of G(y) by a 2 . 
Then Eq. (|lj) gives 

(l-2bt a ) 2 it 2 DV 



-^surv(^) ~" t 



dl exp 



2<tH 2 



I 2 



When t —t 00 the integrand becomes sharply peaked and 
may be evaluated using the saddle point method. This 
gives a — [2z + l)/4 and 

x(6z-l)/4 -Ct (1 ~ 2i)/2 



Psurv (t) ^° t 



(16) 
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FIG. 2: a) Comparison between the survival (filled symbols) and persistence (open symbols) probabilities, b) The corresponding 
local exponents. The horizontal lines correspond to the analytic values given by 9 — 2/(2 — 7). The data for RW survival are 
averaged over variable number of realizations ranging from 10 9 for 7 = to 2 x 10 7 for 7 = 0.5. The DLCA simulations are 
averaged over 50000 simulations on a system of size 55555. The initial distance between particles is 10 (upper curves in Fig. 
a) or 2 (lower curves). 




10" 



FIG. 3: Survival probabilities for 7 = —2 with lo = 
2 (dashed), 3 (dotted), and 4 (dot-dashed). The solid line 
is a guide to eye with a slope /3rw = 0.25. The inset shows 
how the local stretching exponents converge to the same value 
independent of the initial distance lo- 



Inserting the value of a coming from the Lifshitz ar- 
gument to the result of an algebraically expanding in- 
terval, /3dct = 1 — 2a, leads to the same streching ex- 
ponent = (1 — 2z)/2. These two results coincide, 
as a consequence of the peculiar scaling [Eq. (|l5|)] and 
that the tail of the interval length distribution decays as 
G(y) ~ exp(— y 2 ). We emphasize that the fluctuations 
of the surrounding, slow particles determine the stretch- 
ing exponent and that it is purely a coincidence that the 
Lifshitz tail argument gives the same result as the use of 
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t= 10000 




FIG. 4: The scaling plot of the location distribution of the 
rightmost particle in the random walk simulations for 7 = — 2. 
The values of the scaling exponents are z — 1/4 and a = 3/8. 
The solid line shows a Gaussian fit to the data. 



the average value. 

The stretching exponent /3rw = (1 — 2z)/2 has an ob- 
vious interpretation. There are two length scales in the 
problem. The first one is related to the random walkers 
with time-dependent diffusion coefficients, L\ ~ t z , and 
the other to the surviving particle, L2 ~ t 1 / 2 . The argu- 
ment of the exponential decay is simply the ratio of these 
two scales in the problem, P surv (£) ~ exp(— L 2 /Li). Al- 
though reasonable, the calculation above shows the deli- 
cacy of the survival: the distance between the particles 1 
and 3 involves a third, non-trivial length scale, L 3 ~ t a , 
with a — (2z + l)/4. The above considerations can also 
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FIG. 5: The survival probability for 7 < 0. The inset shows 
the bounds for the stretching exponents for the survival (filled 
symbols) and persistence (open symbols). For details see text. 
The dashed [solid] line is given by (1 - 2*0/2 [2(1 - 2z)/3]. 

be made by resorting to an argument which considers the 
two characteristic time scales T\ ~ i 1+7Z and T 2 ~ t. It 
is easy to see, that the ratios between t he scal es obey a 
diffusive like scaling relation L 2 /Li ~ \JT 2 /Ti such that 
any quantity involving the ratio of length scales may be 
given in terms of the ratio of the time scales and vice 
versa. 

In figure |B| the survival probabilities are plotted for 
7 < (for a similar figure for the persistence see fig- 
ure 3 in j[o|). In spite of being able to simulate rather 
small probabilities the asymptotic regime is not reached 
in the simulations. Similar problems with a slow con- 
vergence to the asymptotic value have been encountered 
in other reaction-diffusion systems ]3^, Q and they 
might be overcome by a more efficient use of the cloning 
method (2^, ^6|. The inset of Fig. || shows bounds for the 
stretching exponents as a function of the dynamic expo- 
nent z = 1/(2 — 7). The upper bounds are obtained by 
fitting a line to the three or four largest time points and 
measuring the slope. To obtain the lower bound, we con- 
sidered the change of the local slope and extrapolated to 
1/t — > 0, when it was possible. This method neglects the 
saturation of the local exponent after a finite crossover 
time and therefore gives a lower bound. For compari- 
son, the corresponding bounds for persistence are also 
shown in the inset. There is a clear difference between 
the two. The numerics is consistent with the prediction 
/3rw = (1 — 2z)/2, and for the persistence the data sug- 
gest an expression (5c — 2(1 — 2z)/3. 

The difference between the mean-field model and the 
DLCA is further elucidated in figure ||. It shows that in 
the DLCA the distance distribution between the clusters 
surrounding a persistent one scales similar to that of the 
cluster size distribution 

P(l;t)=L-'h(y). (17) 



FIG. 6: The scaling of the distance distribution between the 
clusters surrounding a persistent one in DLCA for 7 = —0.8. 

Hence, the distribution widens at the same rate as the 
average distance L(t) ~ t z grows in contrast to the RW 
case. For large x the scaling function h(x) ~ exp(— bx) 
and the Lifshitz tail argument leads to an estimate /3l = 
(1 — 2z)/3, which disagrees with the numerics. 

The inconsistency is not surprising since in the DLCA 
there are fluctuations coming from the statistical nature 
of collisions, which are not taken into account in the Lif- 
shitz approach. More precisely, the diffusion constants of 
the neighbors of persistent clusters have some unknown 
distribution. Furthermore, the diffusion constant also 
correlates with the distance from the persistent cluster. 
These facts together with the fact that the stretching 
exponent is determined by the fluctuations makes an an- 
alytical estimation of the persistence for 7 < hard. 

IV. IMPLICATIONS FOR THE CLUSTER SIZE 
DISTRIBUTION 

We now turn to the relation of the persistence to the 
cluster size distribution. We concentrate first on the case 
7 > 0, when the cluster size distribution has a power- 
law tail at small cluster sizes. The dynamical scaling to- 
gether with the definition of the dynamical and polydis- 
persity exponents, z and r, respectively, were discussed 
in the Introduction [Eq. ([]])]. The scaling theory further 
states that all the cluster number densities decay in a 
similar manner at large times, i.e., n s (t) /n\(t) — > b s as 
t — > 00 J35| , where b s is a constant. Here the exponent 
of interest is the universal decay exponent, w, which de- 
scribes the decrease n s (t) ~ t~ w . 

Using S(t) ~ t z together with f(x) ~ x~ T as x — > in 
equation (|l|) gives n s (t) ~ i~( 2 — T ^ z s~ T so that the three 
exponents defined above are related by the scaling rela- 
tion w = (2 — t)z |3(|. Therefore the full characterization 
of the dynamic scaling requires the knowledge of only two 
of the exponents. However, even on the mean-field level 
of Smoluchowski's rate equation theory the only readily 
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TABLE I: Exponents measured from the numerical data. For 
7 = 0.40 the asymptotic regime is not reached in simulations 
(except for z) and only upper bounds are shown. 



7 


z 


e c 


T 


0.00 


0.500 ± 0.001 


1.50 ±0.02 


1.00 ±0.02 


0.40 


0.625 ± 0.001 


< 1.35 


< 0.10 


0.57 


0.699 ± 0.002 


1.43 ±0.05 


0.02 ±0.05 


1.00 


1.00 ±0.01 


2.00 ±0.02 


0.00 ±0.02 



calculable exponent for DLCA is the dynamic exponent 
z. The difficulty with, for example, the polydispersity 
exponent r arises from the fact that to calculate it re- 
quires the knowledge of the whole scaling function p[ . 
Next we argue how knowing the persistence exponent 9c 
helps to overcome this problem. 

Let us start from the trivial size independent case, 
7 = 0, for which an exact solution of the cluster size dis- 
tribution n s (t) is possible [^7|. The actual form of this 
distribution is not important for our purposes. The point 
is that the decay exponent w = 3/2 for any short-range 
correlated initial distribution n s (0). Also the cluster per- 
sistence exponent is universal M. Hence, by noticing 
that for a monodisperse initial condition, n s (0) = <5i jS , 
the persistence probability is simply ni(i), we obtain the 
persistence exponent 9c {0) = w(0) = 3/2. 

The exponents 9c and t should be related also for 
7/0, since the persistent clusters are those ones, which 
have not aggregated with other ones. Asymptotically, 
the number of these clusters will be presented by the 
part s <C S(t) of the cluster size distribution, which in 
turn is characterized by the exponent r. Thus the same 
identification 9c — w can be made also for < 7 < 2 
and we are led to the scaling relation 

9 c = (2-t)z. (18) 

The same relation is valid in a different context of the 
scaling of intervals between persistent regions in the 
reaction-diffusion model A + A — > j3^]. Here 9c = 2z 
and Eq. ( [l8| ) gives r(7) = 0. This is interesting in two 
respects. First, the polydispersity exponent is discontin- 
uous as 7 — > since r(0) = —1^0 = t(0 + ). Although 
quite uncommon, such an outcome is possible also on 
the mean- field level of the rate equation theory |38|]. It 
is more surprising, that the polydispersity exponent is a 
constant, independent of the value of 7. It indicates that 
for any 7 > the physics of small clusters is dictated by 
the fact that they are essentially immobile compared to 
the larger (average-sized) ones in the system. 

Simulations confirm the constant value of t although 
again the crossover effects make the analysis intractable 
near 7 = JTc[ ] - The numerically estimated values for 
the exponents are presented in table |. The scaling rela- 
tion (|l^) is obeyed within the error bars. 

For 7 < the scaling function of the cluster size distri- 
bution behaves as f(x) ~ exp(— Ax** 1 ') when x — » 0. 



Using a similar reasoning as for 7 > leads now to 
the relation (3c = \lA z - Together with the result (3c — 
2(1 — 2z)/3 this suggests that ^(7) = 27/3. Direct mea- 
surement of the exponent fj, is hard as one would need 
to compute the scaling function f(x) for x < 0.1 to see 
the asymptotic behavior. However, even rough numer- 
ics shows that fi(— 2) > —1.75, which is larger than the 
mean-field value predicted by the Smoluchowski's rate 
equation theory, fi = 7. Hence, the spatial fluctuations 
help clusters to survive longer. 

V. CONCLUSIONS 

We have investigated the probability of a cluster to re- 
main unaggregated in one-dimensional DLCA. The diffu- 
sivity of clusters is taken to vary with size as D(s) ~ s 7 , 
which extends the results known for 7 = to the more 
relevant case of size dependent diffusion. 

The first main result is that the persistence probability 
decays as 

!exp(-C^ c ) ,7<0 
t- 3 / 2 ,7 = (19) 

i- 2 /( 2 ^) , 7 >0. 

The stretching exponent fits well to the expression (3c = 
2(1 — 2z)/3 where the dynamic exponent is given by 
z = 1/(2 — 7). Equation ( |l9[ ) shows that one can not 
use the exactly solvable size independent aggregation as 
a starting point for a perturbative analysis of the size 
dependent case. The second main result is that the de- 
cay of persistence is related to the dynamic exponent z 
through the scaling relations 

9 C = (2-t)z 
Pc = \n\z, 

where the exponents r and fi characterize the small size 
tail of the cluster size distribution. Hence, by solving 
for the persistence one determines the behavior of the 
cluster size distribution. For 7 > the scaling relation 
and Eq. (|l9|) lead to a discontinuity of the polydispersity 
exponent: r(0) = — 1 but for < 7 < 2 the distribution 
is flat and t = 0. 

The persistence probability for 7 > is obtained 
from a mean-field analysis for three annihilating ran- 
dom walkers. It explains the discontinuous and non- 
monotonic behavior of the persistence exponent, i.e., why 
3/2 = O (O) > 9 C {0+) = 1. This is since for 7 > 
a persistent cluster eventually adopts the optimal strat- 
egy J28| by becoming more and more stationary as time 
goes on. This interpretation is further supported by the 
fact that the probability of an originally empty site to 
be never occupied by a cluster decays algebraically with 
the same exponent as the cluster persistence ||. The 
major consequence of the discontinuity is the divergence 
of the crossover time to the asymptotic behavior when 
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7 — > + . This also plagues the scaling of the cluster size 
distribution since these two are interconnected. 

The mean-field random walk analysis, which can be 
analyzed in the asymptotic limit when 7 > 0, becomes 
intractable for 7 < 0. We have thus resorted to nu- 
merical studies. These reveal that while the RW picture 
adequately describes the persistence for 7 > 0, it is inad- 
equate for 7 < 0. The reason is that there the persistence 
is affected by the fluctuations in the motion of the slowly 
moving particles around the persistent one. These are 
taken into account approximately in the mean-field the- 
ory, which results only to a qualitative understanding of 
the persistence. For 7 > the approximation is prac- 
ticable as the fluctuations of the slow particle/clusters 
become asymptotically irrelevant. For 7 < they are 
significant as the persistence decays much faster than a 
power law. As an interesting consequence, the mean-field 
theory is applicable when the cluster size distribution is 
broad around the mean (7 > 0: f(x) ~ x~ T , x — > 0) but 
not when it is narrow (7 < 0; f(x) decays faster than 
any power for x — ► 0). 

The difference between the mean-field random walk 
model and the DLCA is demonstrated by the scaling of 
the distribution measuring the distance between the par- 
ticles (clusters) enclosing a surviving (persistent) one [see 
Eqs. ( P^|) and (1171)1 . The main difference is that in the 
theory the average distance grows faster than the dis- 
tribution widens whereas in the DLCA these both take 
place at the same rate. This implies the existence of a 
new, non-trivial length scale ~ t a in the RW-problem. 
A Lifshitz tail argument suggest an expression a = 
(2z + l)/4. This leads to /3 RW = (1 - 2z)/2, which agrees 
with the simulations. Hence, the argument of the expo- 
nential decay is the ratio of the two natural length scales, 
i 1 / 2 and t z , of the problem. An intriguing detail of the 
random walk model is that according to the numerics the 
position distribution of the neighbor of the surviving par- 
ticle scales as p(x;t) = t~ z g ([x — bt a ]/t z ) with a purely 
Gaussian scaling function g(y). It would be worthwhile 
to try to show this analytically and also solve Eq. (]?]) 
with appropriate boundary conditions. This would re- 
quire new analytic tools to handle time-dependent ab- 
sorbing boundary value problems as the traditional im- 
age method can not be applied. We believe this to be an 
unsolved mathematical problem waiting for solution. 



The present study investigates cluster persistence in 
diffusion-limited cluster-cluster aggregation. It would be 
interesting to consider the behavior of unaggregated clus- 
ters in other models, too. Furthermore, we have concen- 
trated only on the one-dimensional case. It is natural 
to ask what can be done in higher dimensions. There 
a similar simple random walk analysis is hardly possi- 
ble. On the other hand the long crossover effects near 
7 = presumably persist and make simulation studies 
hard. Nevertheless, we believe that the general structure 
of the problem remains and conclude with the conjecture 
that also in higher dimensions the behavior of the clus- 
ter size distribution is determined by the solution of the 
cluster persistence problem. 
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APPENDIX A: ASYMPTOTIC ANALYSIS OF 
SURVIVAL FOR 7 > 



The Fourier transform of Eq. <M) reads 



dp 



= -(D 2 + D^ z )(ki + k 2 )p + 2D 2 k x k y p, (Al) 

where we have for notational simplicity used variables 
x and y instead of x\ 2 and X23, respectively. The hat 
denotes the Fourier transform and k x and k y are the 
associated Fourier variables of x and y. The solution 
of Eq. (Al) fulfilling the initial condition pf(x,y;0) — 
5(x - x )S(y - j/o) is 

Pf{k x ,k y ;t) = e ik ^ x e+ lk yV0-D 2 t{k ![: -k y f~V{t)t{kl + k 2 y )^ 



where T>(t) — D\l lz /{p(Z + 1). The subscript / refers to 
the solution without absorbing boundaries. The inverse 
transform reduces to calculating Gaussian integrals with 
the result 



J 



p f (x,y;t\x ,y ) 



*Kty/V(t) [2D 2 + V(t)} 



exp 



[D 2 +V{t)\ [{x - xq) 2 + (y - y ) 2 ] + 2D 2 (x - x )(y - yp) 
4W{t) [D 2 + V(t)] 



(A2) 



At the long time limit this reduces to a Gaussian 

1 ( (x - xq) 2 + (y - y Q ) 2 



pf(x,y;t\x ,yo) 



iwT>(t)t 



exp 



which is nothing but the solution of Eq. ([?]) for D 2 = 0. 
This validates the approximation made in section III C . 
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Since the solution ( [A2] ) is not symmetric in reflection 
with respect to the x- and y-axis, the method of images 
frequently used in problems including absorbing bound- 
aries can not be applied to construct the solution which 
would be zero along the axes. To obtain an estimate for 
the survival probability as a series expansion in powers 
of t, we neglect the cros s-term 2D 2 (x — xo)(y — yo) m the 
exponential of Eq. ( A2) and denote the resulting radially 
symmetric part by pj. The term omitted is of the same 
order in t as the term D 2 [(x — xo) 2 + (y — Do) 2 ] and would 
hence contribute only on the prefactors in the expansion. 
Now the image method gives the solution obeying p = 
along x — and y = for x > and y > 0: 

p(x,y;t\x ,y ) » p s f (x, y; t\x , y ) 

- P S f(x,y;t\ -20,2/0) 

- P S f(x,y;t\x ,-y ) 

+ p S f( x ,y;t\ -x ,-y Q ). 

Integrating this over the first quadrant {x > 0,y > 0} 
yields 



2zx y 
irD 1 t 2z 



[l - 2zK + 6z 2 K 2 + 0(TZ 3 )] , (A3) 



where 7Z = D 2 j (D\t lz ) denotes the ratio of the diffusion 
coefficients. The asymptotic behavior sets in for 1Z <C 1, 
which indicates the divergence of the crossover time to 
the asymptotic behavior when 7^0. 



APPENDIX B: TWO PARTICLE SURVIVAL 

Consider the survival of two particles, which annihilate 
at contact but otherwise evolve according to 

= &(*) 
= &(*). 

where (&(<)) = and (&(<)&(*')) = 2V l {t)6 l] 5{t - t 1 ). 
Let the diffusion coefficients of the particles to be T>\ (t) — 




L»i^ z and V 2 {t) = D 2 , where z = 1/(2 - 7). The dis- 
tance y(t) = X2 (t) — x\ (t) between the particles obeys the 
Langevin equation 



y(t) = y/D^ +D 2 T(t), 

where (T(t)) = and (T(t)T(t')) = 26 (t - t'). This is of 
the standard form [pi] and can directly be transformed 
to a Fokker-Planck equation 



dp{mt) ={D 2 + D 1 t-»f p{mt) 



Of 



dy 2 



(Bl) 



where p(y; t) is the probability density of finding the two 
particles at distance y a t ti me t. 

The solution of Eq. (Bl) fulfilling the boundary and 
initial conditions p(0;t) — and p(y;0) = 6{y — y ) is 
readily found to be 



e -(v-Vo) 2 /4T _ e -(y+ao) 2 /4T 



where T(t) = D 2 t + Di< 1+7Z /(1 + 72). The survival 
probability 



/>oo 

P S urv{t) = dy p(y;t) = erf 
Jo 



yo 

/W 



whose asymptotic behavior at large t is given by 

-fsurv (^) 



' y (TrD 2 t)- 1 / 2 [l-K/{Az) + ... 
y [TT(D 1 +D 2 )t}- 1 / 2 [l-0(t)\ 



^[nD^^z] 1/2 [l-z/K 



.7 < 
,7 = 
• 7>0, 



where 7Z = D 2 / {Dif z ) illustrating again the divergence 
of the crossover time when I7I — > 0. The survival expo- 
nent ^rw = max{l/2, z}, i.e., it is given by the dynamics 
of the faster particle. The interpretation of the result for 
7 7^ is simple: eventually the time scales separate and 
the slower particle becomes stationary. 
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